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d ■ ABSTRACT 

We present results and source counts at 90/im extracted from the 
' Preliminary Analysis of the European Large Area ISO Survey (ELAIS). 

H I The survey covered about 11.6deg^ of the sky in four main areas and 

was carried out with the PHOT instrument onboard the Infrared Space 
Observatory (ISO). The survey is at least an order of magnitude deeper 
than the IRAS lOO/im survey and is expected to provide constraints on the 
formation and evolution of galaxies. The majority of the detected sources 
are associated with galaxies on optical images. In some cases the optical 
associations are interacting pairs or small groups of galaxies suggesting 
the sample may include a significant fraction of luminous infrared galaxies. 
The source counts extracted from a reliable subset of the detected sources 
are in agreement with strongly evolving models of the starburst galaxy 
population. 

Key words: surveys - galaxies: evolution - galaxies: formation - galax- 
ies:Seyfert - galaxies:starburst - infrared: galaxies 



1 INTRODUCTION 

As is well known, a large fraction of the bolometric 
luminosity of galaxies is emitted at far-infrared wave- 



lengths as a result of reprocessing of starlight by dust. 
The spectra of galaxies observed by IRAS typically 
peak in the 60-100/im range. Observations of galaxies 
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in the far-infrared are therefore essential for deter- 
mining the star formation rate and dust masses of 
galaxies. Interest in dust-enshrouded galaxies has re- 
cently been renewed by the discovery of a far-infrared 
(140-850/im) background, thought to be due to dis- 
crete sources (Puget et al. 1996, Fixsen et al. 1998, 
Hauser et al. 1998, Lagache et al. 1999), and detec- 
tions of high redshift galaxies with SCUBA (Smail 
et al. 1997, Hughes et al. 1998, Barger et al. 1998, 
Bales et al. 1999), which are thought to be primeval 
analogues of the starburst galaxies observed locally. 
Deep surveys in the far-infrared therefore promise to 
be very useful for understanding galaxy formation and 
evolution, as they can recover information missed by 
optical/UV studies (e.g. Steidel et al. 1996, Madau 
et al. 1996). Studies of infrared galaxies detected by 
IRAS suggest that violent episodes of star-formation 
are usually associated with very high dust optical 
depths (Condon et al. 1991, Rowan-Robinson & Efs- 
tathiou 1993, Franceschini et al. 1994) which usually 
exceed unity even at 60/im. Models of dusty tori in 
AGN (Pier & Krolik 1992, Granato & Danese 1994, 
Efstathiou & Rowan- Robinson 1995) also suggest that 
the optical depths of type 2 objects are at least as 
high. Observations of galaxies with ISO seem to have 
detected colder dust than observed by IRAS (Kriigel 
et al. 1998, Bogun et al. 1996). 

The European Large Area ISO Survey (ELAIS; 
for a thorough description of the survey and its goals 
see Oliver et al. 2000, PAPER I), that surveyed about 
12 square degrees of the sky at 15- and 90-//m (and 
smaller areas at Q.Tjjbm and 175/im) in mainly four 
(low cirrus) areas spread over the northern and south- 
ern sky, and other ISO surveys (Kawara et al. 1998, 
Puget et al. 1999, Dole et al. 1999) represent the 
first opportunity since IRAS to study the properties 
and evolution of the far-infrared galaxy populations 
at intermediate redshift. An extensive follow-up pro- 
gramme from the radio to the X-rays, already under 
way, ensures that the ELAIS detected galaxies will be 
the subject of thorough study for the next few years. 

In this paper we describe the pipeline developed 
for the production of the ELAIS Preliminary Analysis 
(EPA) 90pm catalogue and present source counts. An 
accompanying paper by Serjeant et al. (2000; PAPER 
II) presents the 6.7 and 15/im pipeline and counts. 
Oliver et al. (in preparation) present an analysis of 
the multiply-covered ELAIS areas. 

The paper is organized as follows: In section 2 
we describe the data reduction and source extraction 
methods used. In section 3 we discuss tests for as- 
sessing the reliability and completeness of the survey 
including comparisons with a parallel pipeline to be 
described in a separate paper (Surace et al. in prepar 
ration). In section 4 we present the source counts. In 
section 5 we discuss briefly our results and outline our 
conclusions. 



2 OBSERVATIONS AND REDUCTION 

As described in detail in Paper I, the data consists 
of a number of pointed observations (in raster mode) 
which use the 3x3 array of the ISOPHOT instrument 
(Lemke et al. 1996) onboard ISO (Kessler et al. 1996). 

The area covered by each raster was typically about 
a quarter of a square degree. The integration time at 
each raster position varied but it was typically 20s. 
The initial observing mode did not allow for any re- 
dundancy. Analysis of some of the early rasters in- 
dicated that the efficiency and reliability of the sur- 
vey could be significantly improved with some redun- 
dancy, so for the latter half or so of the survey we 
switched to an 'overlapping mode' so that each part of 
the sky observed thereafter was covered at least twice. 
Some areas of the survey were observed on more than 
one occasions thus providing further redundancy and 
depth. 

The survey was carried out in four main areas 
(three in the northern and one in the southern hemi- 
sphere) and some smaller areas of special scientific 
interest (PAPER I). 

The pipeline we developed for the Preliminary 
Analysis consists of two basic stages. In the first stage 
we process the raw data using the standard Phot In- 
teractive Analysis (PIA) software (Gabriel et al. 1997) 
until we derive a signal per 2s integration ramp. This 
involves discarding of some readouts at the beginning 
of the ramp, and treatment of cosmic rays. The data 
stream includes data recorded while the telescope is 
slewing between pointings (this typically takes 2 or 3 
seconds per pointing). 

The second stage of our analysis, performed with 
purpose-built IDL routines, involves the extraction of 
point sources from the timeline of each of the 9 pix- 
els of the PHOT array. Given the poor understanding 
of the instrument flat fleld at the time it was consid- 
ered premature to try to extract sources from maps. 
Our source extraction method is therefore very dif- 
ferent from that employed in deep 175/im surveys 
where multiple redundancy and detector stability al- 
low more traditional and reliable techniques to be em- 
ployed (Kawara et al. 1998, Puget et al. 1999, Dole 
et al. 1999). For each pixel we ran an iterative me- 
dian filter that removed outliers (arising from cosmic 
rays) and determined the fractional deviation from 
the background, p, for each ramp. Pointings for which 
the weighted-average of p over the pointing, p, ex- 
ceeded a certain threshold {iop) were fiagged as po- 
tential point sources. 

The pipeline produced source lists on a raster- 
by-raster basis and typically yielded ~ 0.08 detections 
per pointing per pixel. Unfortunately, visual scanning 
of the time profiles of the potential sources revealed 
that only a small fraction of them are likely to be 
real sources. There are a number of effects that can 
give rise to spurious sources and some of these are 
illustrated in Figure 1. 
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Figure 1 shows the data stream of one pixel for 
about 6 pointings centred on a very clear detection of 
a source (at t « 1050 — 1085s) at different stages of 
the processing. Note the delay in the signal response 
to the change in illumination (about 2-3s) which is 
partly due to the movement of the telescope into po- 
sition and partly due to the hysterisis of the detector. 
The latter is also responsible for the gradual fading of 
the signal at the end of the pointing. The occurence of 
these effects at the beginning and end of the pointing 
leaves little doubt as to the reality of the detection. 

The source profiles are usually more complicated 
than the one shown in Figure 1, due to the occurence 
of cosmic ray events. Examples of such events can be 
seen on the same stream of data. Note for example 
the dips in the signal after the spikes at t w 985s 
and 1115s, and the fading profile after a spike at 
t fs 1125s. In addition, the ClOO detectors sometimes 
display a drift behaviour which can give the appear- 
ance of a source if the peak of the cycle coincides with 
the centre of a pointing. 

While a procedure for automatic filtering of false 
detections was developed it was found to be unreli- 
able. The time profiles of all potential sources were 
therefore visually scanned by at least two observers 
in order to construct the final source list. Five differ- 
ent classifications were used ranging from 1 (definite 
detection; see Figure 1) to 5 (severely affected by cos- 
mic rays). 

2.1 Flux Calibration 

The standard method of flux calibration of the PHOT 
instrument is to make use of the internal Fine Calibra- 
tion Source (FCS; Schulz et al. 1998). The FCS has 
been calibrated using planets and asteroids at bright 
fluxes and stars at faint fluxes. The FCS is observed 
immediately before and after each raster in order to 
monitor the change in the rcsponsivity of the detector 
over the course of the raster. Despite the best efforts 
of the instrument team, the calibration obtained in 
this way was considered until recently to be of only 
qualitative value. It is now thought to be uncertain by 
about 20-30% depending on the type of observation. 
This motivated us to pursue an alternative method 
for the flux calibration for the EPA which makes use 
of celestial standards such as the lOO/im background 
measured by COBE/IRAS and IRAS sources. The ap- 
plication of a method based on celestial standards on 
a large and homogeneous dataset such as ELAIS can 
also provide an independent check on the FCS method 
and provide insight into the characteristics of the in- 
strument. 

The basic steps followed in our method can be 
summarized as follows: 

(i) determine the flux per pixel F^^d^t (in Jy) of 
each detection by multiplying p by the predicted back- 
ground flux incident on the pixel F^, = 7^f2, where lu 



is the predicted background intensity and Q is the 
solid angle (0.44468 x lO^^sr; Klaas ct al 1994) sub- 
tended by a ClOO (43.5" x 43.5") pixel. 

(ii) find an empirical conversion factor fp^t for 
converting F^.det to a total flux for a point source 
centred on the pixel Fu,arc- This factor is determined 
by comparing the corrected fluxes with those of cal- 
ibration stars and IRAS sources and includes correc- 
tions for the psf, straylight etc. We will find that fp^t 
is the ratio of two factors fbakg and fpsf- The factor 
fbckg corrects our background estimate for the effec- 
tive solid angle, straylight etc. and fpsf corrects for 
the fact that only a fraction of the flux from a point 
source will be recorded by a single-pixel detection. 

The great advantage of the method is its trans- 
parency. Its validity ultimately depends on the ac- 
curacy with which the background can be estimated 
(given that it is dominated by the zodiacal emission 
which varies with the time of observation) and on the 
accuracy of the fluxes of the sources used as calibra- 
tors. 

In the rest of this section we describe how we 
estimate the background, fbckg and /ps/, and test our 
calibration with stars and IRAS sources. Readers who 
are not interested in the details of these estimates and 
tests can find a summary of our results in section 2.5. 

2.2 Background estimate 

Initially, wo used the COBE background (obtained 
from the IRSKY service at IPAC) which gave us 
very poor spatial resolution. Our approach was subse- 
quently refined to consider the separate contributions 
from the local (zodiacal) background, which depends 
on the time of observation, and Galactic & extragalac- 
tic components. 

The background is estimated by combining the 
contributions of Galactic omission from the lOO/im 
maps of Sclilegel et al. (1998) and a model for the 
zodiacal light. 

The zodiacal light emission depends on the eclip- 
tic latitude fi and the solar elongation e (e.g. CAM 
manual). It is worth noting that due to observing con- 
straints for the satellite the solar elongation of most 
ISO observations was in the range 60-120°. With the 
exception of two small areas (X3 and X6) \(3\ for the 
ELAIS areas lies between 40 and 75°. 

The zodiacal model used is based on the model 
described in the CAM manual (page 77) which gives 
the zodiacal background at 10.9/jm as a function of (5 
and (more crudely) e. The predicted 10.9//m back- 
ground from this model is about 30% higher than 
that measured by ISOCAM towards the ecliptic plane 
(e = 104° and P = -2.4°; Reach et al. 1996) with 
the quoted uncertainty of ISOCAM calibration at the 
time being 15%. 

To extrapolate to far-infrared wavelengths we as- 
sume that the spectrum of the zodiacal light can be 
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Figure 1. Timeline of a typical source with qual=l (definite detection). The top panel shows the raw data (in the form 
of the difference dV between successive readouts in volts) for about 100s centred on the expected time of pointing at the 
source (dotted vertical line). The dashed vertical lines show the time the telescope starts slewing onto a new position. The 
middle panel shows the data at the end of the PIA stage of the data reduction. The bottom panel shows the fractional 
deviation from the background after the median filtering. 



modelled as a blackbody at a (constant) temperature 
of 275K (Hauser et al. 1984). This model gives a ratio 
of lOO/im/25/im backgrounds of 0.16, in very good 
agreement with the estimate of Schlegel et al. (1998). 

In general, the temperature of the dust in the 
zodiacal cloud could vary with /3 and e (as our line 
of sight through the zodiacal cloud samples dust at 
varying distances from the sun) but this is unlikely to 
affect the estimated background by more than 5-10% 
(Schlegel et al.). The range of |/3| and e spanned by 
the ELAIS survey observations is also quite narrow. 

In Figure 2 we compare our background estimate 
, with the flux per pixel determined by PIA^ using 
the FCS calibration, F!^"^. We find that the ratio of 
the two raster-averaged estimates fbckg =< F^'^ > 
/ < Fv > has a median value of 1.7. The origin of 
this difference is unclear but it may be related to the 
incidence of stray light on the detector. What is very 
encouraging is that if we multiply our estimate of the 



* We refer here to the quantity phtiaap.mnfl produced by 
PIA Version 7.3 



background by fbckg, it agrees with that determined 
by PIA within 15%. 

2.3 ELAIS calibration stars 

Given that there are no bright [S > O.GJy) 12/im 
sources in the ELAIS fields (to avoid saturating ISO- 
CAM) we don't expect any 'photospheric' stars to be 
detected at QOfim as these would have a QOfim flux 
< lOmJy . It is possible that there may be some 
stars with circumstellar dust shells, but the useful- 
ness of such stars as calibrators would be limited. In 
order to better determine the ELAIS calibration (as 
well as the general ISOPHOT calibration) three stars 
(HR6132, HR6464 and HR5981) close to the ELAIS 
fields were observed in mini-raster mode (a 3 x 3 raster 
with the star positioned at the centre of a different 
pixel in each pointing). These stars also formed part of 
the ISO ground based preparatory programme. Mod- 
els for their far-IR spectra were constructed by fitting 
near-IR data and extrapolating to longer wavelengths 
as (Hammersley et al. 1998). A more empirical 
approach was given by Cohen et al. (1999). We use 
Hammersley's predictions for the two brighter stars 
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Figure 2. A plot of the ratio of the background obtained using PIA (averaged over all pointings/pixels in a raster) 
< F^'^ > and that from our estimate from a model of the zodia<;al background and maps of galactic emission < > as 
a function of (a) ecliptic latitude /3 and (b) solar elongation. 



but for HR5981 we use the prediction of Cohen as it 
is more rehable for cool stars. The predicted stellar 
fluxes (after convolving with the ISO 90/im filter re- 
sponse) lie in the range (68-311mJy) and so extend 
the flux range to fainter fluxes than possible with the 
IRAS catalogued sources (see section 2.4). The IRAS 
fluxes of these stars agree very well with the model 
up to GO/im. At lOO^im results from SCANPI {alias 
ADDSCAN) for the two brightest stars, which show 
significant detections, also show an excess over the 
models by factors of 2-5. It is not clear whether this 
implies that the stars have an infrared-excess (e.g. due 
to dust) or whether the SCANPI results suffer from 
some artifact introduced in the processing (e.g by cir- 
rus structure). The faintest of the stars (HR5981) was 
observed twice on the same ISO orbit. The integration 
time per pointing in these mini-rasters (40s) is longer 
than that used for the bulk of the ELAIS survey. The 
mini-rasters for the calibration stars were processed 
in the same way as the survey rasters. 

In principle, each detection of each of the stars 
can be used to estimate /p->t, i-e. the correction fac- 
tor for peak flux F^^det to total flux Fu,src for our 
particular observing mode and in the ideal situa- 
tion where a point source is centred in the pixel. We 
find the median value of /p->t to be 2.8. In Figure 



3 we plot the ratio of fp^tF^^det over the predicted 
stellar fiux for each detection. As expected, the em- 
pirically estimated /p^t is ~ ftckg/ fpsf (0,0) where 
fpsf{0, 0) — 0.6 is the theoretically predicted fraction 
of total flux in the central pixel (Klaas et al. 1994, 
Laureijs 1999). 

2.4 Comparison with IRAS sources 

While the ELAIS fields were chosen to avoid strong 
infrared sources, there are a number of IRAS 100/im 
sources in the PSC and FSC which were detected 
in the survey, some of them in a number of pixels. 
The fiuxes of these sources lie in the range 0.5 J y < 
'5(100) < 3Jy and all of them also have 60yum detec- 
tions. Application of the correction factor fp^t, de- 
rived from the calibration observations, to the peak 
fiux generally gave fiuxes too low compared with the 
IRAS fluxes (after correcting for the small change in 
effective central wavelength of the IRAS and ISO fil- 
ters) by factors of 2-3. There are a number of possible 
reasons for this discrepancy: 

(i) Firstly, the optical identifications of these 
sources appear to be galaxy pairs or galaxies with 
sizes comparable to or larger than the pixel size 
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Figure 3. Comparison of ELAIS fluxes witli IRAS fluxes (squares) and tfieoretically predicted fluxes for tlie ealibration 
stars (stars). Tiie 90/im fluxes of tiie IRAS sources are estimated by linearly interpolating between the 60 and 100 firn 
fluxes. The ratio of 90 to 100/im flux is typically 0.8. For the three calibration stars, we plot all the detections (except 
those in the first and last pointings which make the ba<;kground estimation rather diflacult). The scatter in the ratio gives 
a measure of the reproducibility of the results. The predicted fiuxes of the stars are obtained by convolving the models of 
Hammersley et al. (1998), for HR6132, HR6464, and Cohen et al. (1999), for the fainter star HR5981, with the ISO filter 
response function. 



(43.5"). The ISO telescope psf (for the 105/im filter) 
has been approximately fitted with a gaussian with 
FWHM of 50" ± 5" (Heraudeau et al. 1997, MPIA 
internal report). Clearly for such extended sources it 
is not appropriate to scale the peak flux by fp^t- 

(ii) As it was demonstrated by Hacking & Houck 
(1987) in the Deep IRAS survey in the North Ecliptic 
Pole, the 100/ttm fluxes of the survey proper become 
unreliable as the flux limit of the survey is approached 
(^ 0.5Jy at 100/um). The fluxes tend to be overesti- 
mated by factors of up to 2 as they are boosted by 



positive noise fluctuations which are superimposed on 
the true source flux. 

(iii) It is very unlikely that the source (even if it 
were pointliko) would be centred in the pixel having 
the peak flux. In these situations fp^t should be con- 
sidered as a lower limit to the scale factor. 

(iv) Due to the relatively short integration time per 
pointing it is possible that for bright sources the signal 
does not reach its stabilized value, so the flux will be 
underestimated. 
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In order to determine whether a source is ex- 
tended we need to know what fraction of the fiux from 
a point source would fall in an adjacent, or diagonally 
placed, pixel. These fractions were estimated by Lau- 
roijs (1999) for a point source centred in a pixel. His 
results show that fpsf (0, 46) /fpsf (0, 0) = 0.06 and 
/p,/ (46,46)//,,^/ (0,0) = 0.02 where fj,sf{x,y) is the 
fraction of the flux in the pixel centred at x,y arc- 
seconds from the point source. This assumes all the 
pixels have the same roll auglo. Note that although 
each ClOO pixel is 43.5" square, the separation of pixel 
centers is 46" (Klaas et al. 1994). 

According to this analysis we can assume that 
adjacent detections with flux greater than ~ 15% of 
the peak flux correspond to extended sources or mem- 
bers of group of sources. There are, of course, situ- 
ations where this criterion will fail, that is when a 
point source lies at the corner or edge of a pixel. We 
have simulated this situation and found that we are 
unlikely to overestimate the flux by more than 30%. 
Laureijs also calculated the fraction of the psf in the 
whole array and found it to be ^ 0.8 for the C90 fll- 
tcr, which also sots a limit of ~ 30% to the degree 
by which we would overestimate the flux of a point 
source in this case. 

2.5 Summary of calibration method 

In summary, our strategy for flux calibration is as fol- 
lows: for sources with more than 1 reliable detections, 
the total flux is calculated by j '^^(o'o) X^i -^Arfet> 
where F^^^^ arc the fluxes of non-overlapping detec- 
tions above 15% of the peak. The factor fbckg is esti- 
mated to be 1.7. The factor fpsf{0, 0) is assumed to be 
0.6 (Klaas et al 1994). For single detections we calcu- 
late the flux as fbckgF„^det/ < fpsf > where < fpsf > 
is estimated from simulations of a point source with 
a FWHM of 50" at different positions in the pixel to 
be ~ 0.4. 

In Figure 3 we compare the predicted 90fj,m 
fluxes of the IRAS sources with the fluxes calibrated 
in this way. In conclusion, our approach to the flux 
calibration gives fluxes that agree with those of IRAS 
sources and bright stars (which bracket the flux range 
covered by ELAIS) within 30-40%. 

3 RELIABILITY AND COMPLETENESS 

3.1 Comparison with an independent 
pipeline 

We have pursued a programme of detailed compar- 
isons of the IC detection lists with those compiled by 
a pipeline developed in parallel (but somewhat later, 
to take advantage of a better understanding of the 
instrument) by the group at the Max-Planck-Institut 
fiir Astronomic (MPIA) in Heidelberg. Results from 
that pipeline will be described in detail by Surace et 



al. (in preparation). The detections obtained by the 
MPIA pipeline were also eyeballed with the same clas- 
sification scheme applied to the IC detections. 

The sample of sources on which the preliminary 

90/irre counts are based is a reliable subset of the EPA 
list which was constructed as follows: 

(i) A list of detections extracted by either pipeline 
and which received the top two eyeballing classifica- 
tions (1 or 2) by at least 2 observers was first con- 
structed. 

(ii) The combined list was then trimmed to include 
only detections which received one of the top two clas- 
sifications by at least 3 observers (275 detections). 

(iii) The detections were then merged with a near 
neighbour algorithm (search radius of 1') to produce 
a source list consisting of 163 sources. Only one of 
the sources with S > lOOmJy was detected by the 
MPIA pipeline only; otherwise, the lists from the two 
pipelines are in complete agreement above this limit. 

In Figure 4 we compare the fluxes of the QLA 
reliable detections {not sources) as determined by the 
method described here (i.e. fbckgF^^det) with those de- 
rived from the MPIA pipeline. While there is consid- 
erable scatter (comparable to the scatter seen in the 
comparison of the backgrounds) it is rather encour- 
aging that above 50m Jy the fluxes agree within 40%. 
At the bright end the MPIA fluxes may be under- 
estimated because of self-subtraction. Below 50mJy 
there seems to be a group of detections where the 
fluxes derived by the two methods differ by factors of 
2-3. 

3.2 Simulations 

The data are strongly affected by cosmic rays and 
other detector effects so the coverage is very inhomo- 
geneous. In order to extract the source counts we need 
to quantify the incompleteness as a function of flux 
and position in the survey area. There are two main 
causes of incompleteness. Firstly, sources that lie at 
the edges or corners of pixels are likely to be missed. 
In areas where there is redundancy (either because 
the overlapping mode is used or because the area has 
been re-observed) the probability of missing sources 
decreases. Secondly, some genuine sources may have 
been rejected at the eyeballing stage because they co- 
incide with cosmic rays etc. 

To make a rough estimate of the degree of in- 
completeness due to these (and possibly other) effects 
we have simulated a number of rasters and processed 
them in an identical fashion as the survey rasters. 
The simulated data are constructed by adding scans 
of randomly placed sources onto real data extracted 
from several rasters (where the rasters and pixels are 
chosen at random). The detected sources were merged 
in the same way as the real sources. For practical rea- 
sons the simulated detections were eyeballed by a sin- 
gle observer (AE). This could introduce a bias in the 
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Figure 4. Compajrison of the fluxes of the QLA reliable detections with those derived from the MPIA pipeline. 



completeness correction as we require confirmation 
(i.e. eyeballing class 1 or 2) by at least 3 observers for 
a detection to appear in our final list of real sources. 
To quantify the effect of this we can use the eyeballing 
experience of the whole survey to find the rmrnber of 
detections confirmed by AE which finally appeared in 
the combined list (186 out of 286). For comparison, 
the number of detections not confirmed by AE which 
appeared in the combined list is 34. As the number 
of detections confirmed by AE only (286), and by at 
least 3 observers (186+34), are comparable we con- 
clude that there is not a significant bias introduced 
by using the eyeballing results of a single observer for 
the purpose of estimating the completeness. In order 
to simulate regions of multiple redundancy we simu- 
lated rasters with different steps in one of the raster 
dimensions. Two possibilities were simulated: half an 
array step, to simulate regions of single redundancy, 
and quarter array step, to simulate regions of triple 
redundancy. The sources are calibrated in the same 
way as the real sources. 

The results of this study are plotted in Figure 
5. The fraction of sources detected increases signif- 
icantly as we move from no redundancy to single 
redundancy. Further redundancy improves the com- 
pleteness at faint fluxes. 



4 SOURCE COUNTS 



4.1 IRAS counts 

We have extracted 90/irre cxtragalactic counts at 
fluxes brighter than possible with ELAIS by using 
the PSCz catalogue (Saunders et al. 2000). We ap- 
plied a cut with galactic latitude (|6| > 20) which 
limits the areal coverage to 7.642 storadians (Rowan- 
Robinson et al. 1991) and the number of galaxies to 
11405. We also excluded galaxies with low IRAS flags 
{fqual < 3) at 100/xm. The resulting sample consists 
of 8529 galaxies. The 90/xm flux of each galaxy is es- 
imated by linearly interpolating between the 60/im 
and lOOyum fluxes (by design all galaxies in the PSCz 
catalogue have reliable 60/im fluxes) . 

The IRAS counts are plotted in Figure 7. The 
structure in the counts above 5Jy, which is also evi- 
dent in the 60/im counts (Hacking, Condon & Houck 
1987), is probably due to a local excess in the source 
density dominated by the Virgo cluster. 

The PSCz catalogue only includes galaxies with 
S(60) > O.GJy. IRAS galaxies are known to have 
S'( 100)75(60) ratios at least as high as 3.5 (e.g. 
Rowan- Robinson & Crawford 1989, Efstathiou et al, 
in preparation). A 90/im sample drawn from PSCz 
is therefore bound to be incomplete below 2Jy. This 
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Figure 5. Detection fraction of point sources as a function of flux and redundancy (diamonds=3, stars=l, crosses=0) 
estimated from simulations. The uncertainties are the statistical /(I — f)N errors on the detection fraction / where N 
is the number of sources in each simulation. 



incompleteness is most probably the reason for the 
flattening of the IRAS counts between l-2Jy. 



4.2 ELAIS counts 

To extract the source counts wc need to determine 
the area over which the survey is sensitive to as a 
function of flux. To do that we use the results from 
the simulations given in Figure 5. While we have not 
simulated sources with S > 0.5Jy we can reasonably 
assume that we arc close to 100% complete for those 
fluxes as we have detected all IRAS sources in the 
flelds. We only consider regions of redundancy 3 or 
less in this paper (amounting to a total of lldeg^) 
and therefore exclude the S2 area. All of the small 
areas of special scientific interest arc also excluded 
as some are centred on known objects. The resulting 
areal coverage as a function of flux is given in Fig- 
ure 6. The counts are corrected for incompleteness by 
normalizing each source by the area over which each 
source could be found (Hacking & Houck 1987, Oliver 
et al. 1997). The resulting integral counts are given 
in Figure 7. Flux uncertainties of ±0.15dex are also 
indicated. The uncertainties in the number counts are 
the Poisson errors, uncertainties in the effective area 
have not been included. 



The sources presented here axe extracted from 
a reliable subset of the EPA sources as described in 
section 3. About 42% of the sources brighter than 
140mJy arc detected in at least two neighbouring 
or overlapping pixels. About 67% of the sources also 
have an association with a mid-infrared ELAIS source 
within a 1 arcminute radius. Given the relatively low 
source density of ELAIS sources the probability of 
a chance association is only a few per cent. We also 
note that there arc good optical ids associated with 
the majority of the PHOT sources. 

To obtain an estimate of the likely contamination 
of our sample by cirrus we use the models of Gautier 
et al (1992). Our observing mode is different from 
either of the two cases considered in that paper as 
our reference (background) position is offset by 67" 
to 390". If we use their double aperture mode with 
this range of offset and assume the slope of the power 
spectrum of the cirrus is in the range -2.6 to -3.8, then 
the expected rms cirrus confusion noise is predicted to 
be in the range 0.2-3.1mJy for a background intensity 
of IMJy/sr. The latter is about a factor of 2 brighter 
than that found in the typical ELAIS area (PAPER 
I), so it is very unlikely that any of our sources are 
due to cirrus. 

We estimate the slope of the counts (assuming 
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Figure 6. Area (in square degrees) over whieh the sensitivity of the survey is S. The area is estimated from simulations 
for different values of S and different degrees of redundancy. 



they can be described by a power-law) by applying 
a minimization procedure (the IDL routine linfit) 

on the model 

dN 

In the range 0.158-lJy we find a = 3.71 ± 0.13 

and b — —1.92 ± 0.23. If wo use only sources found in 
areas with redundancy in our analysis, the estimated 
parameters are a = 3.65 ± 0.22 and b = —2.11 ± 0.34. 
We also estimate the slope of the IRAS counts in 
the range 2-20Jy (where the efi^ects of incomplete- 
ness and large scale structure arc minimized) to be 
a = 3.73 ± 0.02 and b = -1.48 ± 0.03. The ELAIS 
counts therefore show some evidence for departure 
from the Euclidian slope. However, because of the 
limited statistics of the ELAIS survey and the in- 
completeness of IRAS at around IJy it is difficult to 
constrain from this study the flux at which this takes 
place. 



5 DISCUSSION 

5.1 Comparison with evolutionary models 

The recent dramatic improvement in observational 
constraints in the star-formation history of the Uni- 
verse and therefore the evolution of the starburst 
galaxy populations in the optical/UV, as well as re- 
cently the submillimeter, has stimulated the develop- 
ment of evolutionary models (Pei & Fall 1995, Pear- 
son & Rowan-Robinson 1996, Pranceschini et al. 1997, 
Guiderdoni et al. 1998, Rowan- Robinson 1999). In 
Figures 7 and 8 we compare the observed counts with 
the models of Rowan- Robinson, Guiderdoni et al. and 
Pranceschini et al.. The counts predicted by the mod- 
els differ by a factor of 3-4 at about lOOmJy and the 
ELAIS counts can potentially discriminate between 
them. 



The model of Guiderdoni et al. is set within the 

framework of hierarchical growth of structures accord- 
ing to the cold dark matter model, and extends earlier 
studies (e.g. by Cole et al. 1994) to the IR/submm 
wavelength regime. In Figure 7 we plot the prediction 
from their models A and E. The latter model incor- 
porates a heavily extinguished (ULIRG) population, 
assumed to dominate at high redshift in order to ac- 
count for the far-ir background. 

The starting point in the model of Rowan- 
Robinson (1999) is an assumed star formation history 
of the Universe, described by an analytic formula in- 
volving 2 free parameters {P,Q). The star formation 
rate is assumed to evolve due to pure luminosity evo- 
lution. The parameters P and Q are then constrained 
by fits to the observed star formation history, the ob- 
served 60- 175 and 850/im counts and the observed 
far-ir background. The model assumes a number of 
infrared galaxy populations which are modelled with 
the radiative transfer models of Efstathiou, Rowan- 
Robinson & Siebenmorgen (2000). The prediction of 
the Rowan- Robinson model is compared with the ob- 
served counts in Figure 8. 

The model of Pranceschini et al. (in prepara- 
tion) assumes that the extragalactic population is 
composed of two components with different evolution 
properties: (1) a non-evolving galaxy population in 
which the far-IR flux is mostly contributed by evolved 
stars; (2) a population of strongly evolving starbursts, 
with peak cmissivity around redshift 0.8 and roughly 
constant cmissivity at higher z. The evolution rate 
for population (2) is optimized to reproduce the mid- 
IR counts and redshift distributions. The model of 
Pranceschini is also plotted on Figure 8. 

While the observed counts tentatively favour the 
Guiderdoni model, the uncertainties due to flux cali- 
bration and correction for incompleteness do not al- 
low us to definitely discriminate between the mod- 
els at the present stage. A more detailed analysis of 
the ELAIS 90/im survey (already under way) should 
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be able to discriminate between the models. Further 
discrimination of the two models should be possible 
with the redshift distributions of the detected galax- 
ies. The models of Pranceschini et al. and Rowan- 
Robinson predict that the vast majority of sources 
at S > lOOmJy will be at redshifts less than 1. 
The Guiderdoni et al. model overpredicts the number 
of high redshift galaxies in the North Ecliptic Pole 
(NEP) survey, although the counts in the latter are 
known to be dominated by a supcrclustcr (Ashby et 
al. 1996). The ELAIS survey weis specifically designed 
(PAPER I) to minimise the effects of large scale struc- 
ture. 

5.2 Comparison with other Infrared and 
submillimeter surveys 

Until recently the deep IRAS 60/itm in the NEP sur- 
vey (Hacking & Houck 1987) provided the deepest ex- 
tragalactic source counts in the infrared. Augmented 
with counts from shallower surveys (e.g. Rowan- 
Robinson et al. 1991) they provided the basis for stud- 
ies of the evolution of the starburst galaxy population 
(Hacking et al 1987, Oliver et al. 1992, 1995, Pear- 
son & Rowan- Robinson 1996). Unfortunately redshift 
measurements of the redshift sources, which should 
have been able to constrain further the form of evo- 
lution, discovered the presence of a 2: = 0.088 super- 
cluster in the NEP field (Ashby et al. 1996). 

Oliver et al. (1997) confirmed the strong evolu- 
tion seen in IRAS surveys with deep counts at 6.7 
and 15/im in the Hubble Deep Field. Counts at 15/jm 
well in excess of no-evolution expectations were also 
reported by Elbaz et al. (1999). Serjeant et al. (2000, 
Paper II) found the evolving high-z population dom- 
inates over the local counts at 5(15) < lOmJy. 

With the advent of SCUBA, submillimeter sur- 
veys (at 850/xm) which exploited the negative K- 
correction due to the spectral shape of starburst 
galaxies, detected source densities 2-3 orders of mag- 
nitude higher than expected from non-evolving mod- 
els (Small et al. 1997, Hughes et al. 1998, Barger et al. 
1998, Eales et al. 1999). While the process of identify- 
ing and determining the properties of those SCUBA 
sources is still incomplete, it appears that a rmmbor 
of them are similar to ULIRGs at 2 ~ 2 - 4. SCUBA 
surveys appeax to have resolved 94% and 34% of the 
COBE-FIRAS background values at 850 and 450Atm 
respectively (Blain et al. 1999). 

FIRBACK 175/^m counts in the Marano 1 field 
also appear to have shown an upturn in the counts 
(Puget et al. 1999). These counts are 3 times higher 
than the models of Guiderdoni et al. and Franceschini 
et al. However, 175/irre counts in a larger area (Dole et 
al. 1999) suggest that the counts in the Marano field 
may have been anomalously high. It is also not clear 
whether the excess in the counts, if any, is due to high 
redshift staxbursts or local cold galaxies. Counts 3-10 
times higher than expected from no-evolution mod- 



els were also found by Kawara et al. (1998) in their 
175/im survey of the Lockman Hole. More results are 
expected from FIRBACK and ELAIS surveys at this 
wavelength. Brighter counts are also expected from 
the ISO Serendipity survey (Bogun et al. 1996, Stickel 
et al. 1998). 

The extensive multi-wavelength coverage and 
follow-up programme of ELAIS ensures that, unlike 
other far-ir surveys where resolution and sensitivity 
can be a serious obstacle, the 90/im galaxies will be 
thoroughly studied (e.g. Morel et al. in preparation) 
and will therefore provide a firm basis for studying 
the evolution of dust-enshrouded galaxies. 

A more detailed analysis of the ELAIS PHT sur- 
vey is underway and results will be presented in sub- 
sequent papers in this series. 



CONCLUSIONS 

We have presented 90/im source counts extracted 
from the catalogue produced by the ELAIS Prelim- 
inary Analysis and PSCz redshift survey. All but 
one of the sources brighter than lOOmJy have been 
confirmed by an independent pipeline developed at 
MPIA, Heidelberg. We compare our estimated fluxes 
with IRAS fluxes and models for standard stars and 
estimate the uncertainty in flux calibration to be 
about 30-40%. 

The ELAIS counts extend the IRAS counts by an 
order of magnitude in flux. The slope of the counts 
in the 0.158-lJy range shows some evidence for de- 
parture from the Euclidian slope. This is consistent 
with the strong evolution seen in other infrared and 
submillimeter surveys. 

Within the uncertainties Eissociated with the flux 
calibration of the survey, the counts agree with the 
strongly evolving models of Rowan- Robinson (1999), 
Guiderdoni et al. (1998) and Franceschini et al. We 
expect the redshift distributions arising from the 
spectroscopic follow-up of the survey to allow us to 
discriminate between these models. 
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Figure 7. ELAIS/IRAS 90/im source counts. The solid and dotted lines are the counts predicted by the models A and E 
(respectively) of Guiderdoni et al. (1998). 
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Figure 8. ELAIS/IRAS 90 fim source counts. The solid lino gives the counts predicted by the model of Rowan- Robinson 
(1999). The dotted line shows the counts predicted by the model of Franceschini. 
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